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Abstract Multiple solutions exist in various experimental situations whenever the sum 
of several amplitudes is used to fit the experimentally measured distributions, such as the 
cross section, the mass spectrum, or the angular distribution. We show a few examples 
where multiple solutions were found, while only one solution was reported in the publica- 
tions. Since there is no existing rules found in choosing any one of these solutions as the 
physics one, we propose a simple rule which agrees with what have been adopted in previ- 
ous literatures: the solution corresponding to the minimal magnitudes of the amplitudes 
must be the physical solution. We suggest test this rule in the future experiments. 
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1 Introduction 

In quantum mechanics, a physics observable is proportional to the modulus of the 
amplitude squared. In case of more than one amplitude contributing to a process, they 
are summed to obtain the total amplitude (generally there are relative phases between 
these amplitudes), and thus one generally has contribution from interference term to the 
physics observable. It is simple to predict a physics observable when the amplitudes and 
the relative phases between them are known, since there is no ambiguity in this procedure. 

However, in many circumstances, the experimental quantities are measured, and from 
which we extract the information on the amplitudes. As there is a square operation 
between the observable and the amplitudes, we would expect multiple solutions in solving 
the equation from a pure mathematics point of view. Then we face a common problem: 
which solution is the physics one. 

In this Letter, we present a few examples where two-solutions were reported in fitting 
with the coherent sum of two amplitudes; and revisit a few cases where only one solution 
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was reported and we find the other distinctive solution in fitting with two amplitudes; we 
also discuss more complicated situations where more than two amplitudes are used to fit 
the data. Finally, we propose a conjecture on how to choose from multiple solutions the 
physics one. 

2 Examples with two-solutions reported 

A few recent examples reporting multiple solutions are in the study of the so-called 
Y states via initial state radiation (ISR) by the Belle experiment. Figure [T] shows the 
invariant mass distributions of 7t + tt~ J/ip and n + 7i~tfj(2S) after all the selection in Belle 
data pfl[2], together with a fit with two coherent resonant terms and an incoherent back- 
ground term. Table [1] shows the fit results, including the F(4008) and F(4260) from the 
7T + 7T~J/ip mode, and the V(4360) and F(4660) from the 7r + 7r~ip(2S) mode. It should be 
noted that in both channels there are two solutions of exactly the same goodness-of-the-fit, 
with exactly the same mass and width for the resonances but with very different coupling 
to e + e~ pair (r e + e -). Instead of choosing one from the two solutions, Belle reported both 
in their publications, however, Particle Data Group (PDG) quoted sometimes only one 
of the solutions in averaging with those from other experiments [3]. Such a treatment of 
data is rather suspicious, since the solutions picked up randomly from different experiment 
may have distinctive features and should not be averaged together. 

Another example is the study of the decay dynamics of rf — > 77r + 7r~ mode. When 
the 7T + 7r~ invariant mass distribution is fitted with coherent sum of the p resonance and 
a contact term, it is found that there are two solutions of equal goodness-of-the-fit. One 
solution corresponds to constructive interference between the two amplitudes while the 
other destructive interference [4]. Similar study in previous analyses only reported one 
solution [5] which corresponds to the solution with constructive interference in the BES 
data [1]. 

3 Examples with one-solution reported 

As the examples shown in the previous section are basically a fit to an s-dependent 
distribution with two coherent amplitudes, one would expect that there would be in 
general two solutions in such a circumstance. In the literatures, there are many such 
kinds of fits, especially in the low energy e + e~ annihilation experiments, where the cross 
section of e + e~ — > hadrons is parameterized as the coherent sum of the amplitudes of 
the vector mesons. We show two typical examples below where only one solution was 
reported, and we redo the fit to obtain the other solution from the data. 
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Figure 1: The ir + n J/ip (upper) and tt + tt ip(2S) (lower) invariant mass distributions and 
the best fit with two coherent resonances together with a background term. The data are 
from Belle [IJ[2]. 



3.1 Branching fraction of — ► utt 

The most precise data on e + e~ un near the <fi resonance were reported by KLOE [6] 
with both uj — > 7t + 7t~7t° and uj — > 771" . The cross section as a function of yfs is parame- 
terized as 

o-(a/s) = a nr (y/s) ■ 

in the KLOE paper [6], where o" nr (y / i) = cr + a'(y/s — M^) is the bare cross section for 
the non-resonant process, parameterized as a linear function of y/s; Z is the interference 
parameter, while M$, and = — s — iM^Y^ are the mass, the width, and the 
inverse propagator of the <fi meson, respectively. 

We take the KLOE data from Table I of Ref. [6] with u — > 7r + 7r~7r° and fit with the 
same parametrization given above. In our fit, the Born-order cross section is calculated as 
Ovis/<5rad, and only statistical errors are considered in the x 2 construction. Table [2] shows 
the results from our fit, where two solutions are found with the same fit quality. We can see 
that the parameters from Solution I are quite similar to those listed in Table II of Ref. [B] , 
the slight difference is due to the simplified procedure for the Born-order cross section 
calculation in our fit. It is found that Solution II has a much larger resonant amplitude, 

o 



1 -Z 



M rh r„ 



Table 1: Fit results of the 7r + 7r~ J/ip and n + n~ip(2S) invariant mass spectra. The first 
errors are statistical and the second systematic. M, T to t, and B ■ T e + e - are the mass (in 
MeV), total width (in MeV), product of the branching fraction to hadronic mode and 
the e + e~ partial width (in eV), respectively. <fi is the relative phase between the two 
resonances (in degrees). 



Parameters Solution I Solution II 





M(y(4008)) 


4008 ± 40t^ 4 




r tot (F(4008)) 


226 ± 44 ± 87 


B 


•r e+e -(F(4008)) 


5.0 ±1.4t^ 12.4 ± 2.41I 4 ! 8 




M(y(4260)) 


4247 ± 12+H 




r tot (Y(4260)) 


108 ± 19 ± 10 


B 


•IW(F(4260)) 


6.0 ±1.2tol 20.6 ±2.3^;} 







12 ±29^ 8 — 1 1 1 ± 7±§f 




M(F(4360)) 


4361 ± 9 ± 9 




r tot (F(4360)) 


74 ± 15 ± 10 


B 


•r e+e -(F(4360)) 


10.4±1.7±1.5 11.8 ±1.8 ±1.4 




M(F(4660)) 


4664 ± 11 ± 5 




r tot (F(4660)) 


48 ± 15 ±3 


B 


•r e+e -(F(4660)) 


3.0 ±0.9 ±0.3 7.6 ±1.8 ±0.8 
39 ± 30 ± 22 -79 ± 17 ± 20 



and the resulting branching fraction of the isospin violating process — > lott is at per 
mille level, about two orders of magnitude higher than Solution I, namely, the solution 
reported in the original work [6]. Figure [2] shows the fit results and the contribution of 
each component in the fit. Similarly, we try a fit to the cross sections measured with 
u — > 77T°, there are also two solutions obtained from the fit and the results are in good 
agreement with those from uj — > 7r + 7r~7r°. 

3.2 F n and the p-u mixing 

Even before the discovery of the r lepton, Tsai calculated the branching fractions of 
such a heavy lepton decaying into vector final states by the conserved vector current (CVC) 
hypothesis using e + e~ — > tt + tt~ data measured in e + e~ annihilation experiments [7j. This 
calculation, referred to as B^Y^.-^ , was tested ever since the discovery of the r lepton. 
As the precisions of both e + e~ annihilation and r decay experiments improved significantly 



Table 2: Results from fits to the e + e — > ujtt cross sections measured with u — > 7r + 7r 7r°. 



Parameter Solution I Solution II 

(T [nb] 7.88 ± 0.04 7.88 ± 0.08 

3?(Z) 0.106 ±0.004 0.106 ±0.006 

Q(Z) -0.103 ±0.003 -1.90 ±0.006 

a' [nb/MeV] 0.064 ± 0.002 0.064 ± 0.006 

S(0 -> wvr ) 4.61 x 10~ 5 7.62 x 10~ 3 



in the past three decades, the test has reached the precision of better than 1% level, and 
the discrepancy of 2cr level between B^Y^-^o^ and B T -^- n o UT is observed [8]. Many 
theoretical efforts have been made to understand this discrepancy, including a better 
understanding of the isospin violation correction and so on. 

One source of the isospin breaking correction is the p-uo mixing effect in the e + e~ 
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annihilation data which is absent in the r decays. The mixing effect on B^ 
estimated by fitting the e + e _ — > 7r + 7r~ data in the vector meson dominance (VMD) model 
(including p and its excited states and uj) and subtract the uj contribution by setting its 
amplitude to zero. This has been carried out in many previous analyses [8j [H [10] , but 
in all cases only one solution was reported. As a heuristic example, we only use the 
CMD2 data and follow the fit described in the CMD2 paper [TT], namely, with the GS- 
parametrization [12] of the F n : 



where 



^GS, 



BW 



GS/ 



BWf (s) 



rriy(l + d ■ T v /m v ) 



ni 



v 



s + f(s) — im v T v (s) 



m,. 



mt. 



s - im u T u 



The definitions of all the quantities can be found in Refs. p^Ol [12] . Here, complex numbers 
S and j3, as well as the mass and width of the p, are fit parameters. The masses and widthes 
of uj and p' are fixed to their PDG values [3], and the phase of (3 is fixed to 180°. 

We found two solutions in our fit, as shown in Fig. [3] and Table [3j It is clear that 
Solution I is in good agreement with the results reported by CMD2 experiment [IT] , and 
the resulting p-u mixing corrections to B ( ^Y^ n - 7T o UT and a M are in agreement with the 
results in Refs. [8j [9]. However, the relative strength between u and p is different in the 
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Figure 2: Fit to the e + e~ — > ujtt° cross sections as a function of center-of-mass energy. 
Data points are shown with dots with error bars. Top plot is for solution I: at y/s = 
1.019 GeV, from top to bottom, the solid curves are continuum, best fit, resonance, and 
interference terms. Bottom plot is for Solution II: at y/s = 1.019 GeV, from top to 
bottom, the solid curves are resonance, continuum, and best fit, the interference term is 
not shown. 

other solution, and the corrections to - o„ and a,, are different too, as shown in 

Table 

As the amplitude of the oj term in Solution II is more than an order of magnitude 
higher than that in Solution I, we would expect a branching fraction of uo — > 7r + 7r~ much 
larger than that reported in previous analyses: about 70% in Solution II versus 1.5% in 
Solution I. It seems Solution II is in contradiction with our expectation that the isospin 
breaking decay oj — > 7r + 7r _ should be much smaller than the isospin conserving decay 
uj — > 7T + 7r~7r°. However, it is obvious that we are used to the small u — > tt + tt~ branching 
fraction because we have never thought about the possibility of the existence of the second 
solution in the fit. 



4 Examples with more than two amplitudes 

The examples shown above are described by the sum of two amplitudes with a free 
relative strength and relative phase. In the case of more than two amplitudes contribut- 
ing to the distribution, the number of solutions is 2 n_1 , where n is the number of free 
amplitudes. 

In the fit to the e + e~ — > 07r + 7r~ cross sections |13j . incoherent and coherent sums of two 
amplitudes (the 0(1680) and K(2175)) are tested, while in estimating the significance of 




s (GeV^) s (GeV^) 

Figure 3: Fit to the e + e~ — > 7r + 7r~ form factors below s = 1 GeV 2 measured at CMD2 [TO]. 
Top: along the data points (dots with error bars) is the best fit. The 1 = 1 part of the fit 
is shown for the two solutions. Bottom: comparison of the p-uj interference part in two 
solutions. 



a third resonance pf (2400)), Belle also tested the possibility of the coherent sum of three 
amplitudes. This procedure is repeated with combined BaBar [T4] and Belle data [13] in 
Ref. p. 5]. In this latter case, four solutions are found with the same masses and widthes 
for the resonances, but with different coupling constants. 

The same phenomenon was also shown in the fit to the combined BaBar [16] and 
Belle [2] data on e + e~ — > 7i + 7i^ip(2S) [17], where two pairs of solutions were found when 
fitting with the coherent sum of the F(4360), F(4660), and y(4260). 



Table 3: Results from fits to the e + e — > 7i + n form factors measured at CMD2 [TO] . 



Parameter 


Solution I Solution II 


Davier [8] 


m p [MeV] 


775.9 ± 0.5 




T p [MeV] 


146.0 ±0.8 




|*| [xlO- 3 ] 


1.62 ±0.06 21.97 ±0.04 




fa [°] 


10.1 ± 1.4 86.56 ±0.17 




101 


0.086 ± 0.004 




Ai3 mixing [ % ] 
Aa mixing [XQ-iO] 


-0.03 ±0.01 ±0.04 ±0.01 
±2.5 ±0.2 ±1.6 ±0.2 


-0.01 ±0.01 
+2.80 ±0.19 
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5 Partial wave analysis 



All the examples discussed so far are one-dimensional problems. In the case of fitting 
a multi-dimensional distribution, the same multi-solution problem also exists. The partial 
wave analysis (PWA) and the Dalitz plot analysis techniques used widely in hadron physics 
potentially have such a problem. 

In a recent BES analysis of the ip(2S) — ► 7r + 7r~J/?/> process [18], the 7r + 7r~ system is 
fitted with contributions from a, a contact term, and a possible tensor amplitude f2- It 
is curious to us why the two very similar models shown in Figure 4 of Ref. [18] give very 
different fractions of the a term. The aforementioned expositions lead us to suspect that 
two different solutions are found in the two models, one with constructive interference 
and the other with strong destructive interference. The authors may try to excavate the 
other solutions in each of those models used in fitting the experimental data. Fortunately, 
according to the study above, we found that all different solutions result in identical line 
shape of the resonance (the same mass and width), so the resonant parameters should be 
correct even if the branching fractions obtained from these analyses are not reliable. 

A general partial wave analysis may have many amplitudes and all of them are added 
coherently to fit the data distributions. In principle, there are multiple solutions, although 
we cannot judge how many fold of ambiguity. In some circumstances, when more than one 
solutions are found, a simple average is taken to be the best estimation of the parameters, 
this is obviously non-physical, since any one of the solutions is a good description of the 
data but the average is generally not; while in some other cases, this is just overlooked, 
which makes the comparison between experiments looks very strange, since one compares 
two different solutions from two experiments which are different by definition. Herein, 
great care must be taken in treating the magnitudes of the amplitudes in extracting them 
from multi-dimensional distributions with the coherent sum of many amplitudes. 

6 Minimal amplitude conjecture 

In preceding sections, we examined a few cases where amplitudes are extracted from 
experimentally observed distributions, and in all the cases more than one solutions were 
found. We also discussed the similar problem in PWA. All these indicate that this problem 
actually is common in data analysis. Up to now, this problem has not been taken into 
consideration seriously. Moreover, no one has thought about how to pick up one solution 
from all the possibilities, and there is indeed lack of solid physics argument on how to 
choose one of the solutions as the physics one. 

However, from the existing examples and our observation, we notice that the solution 
with the smallest modulus of the amplitude agrees with present expectation, such as in 
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the — > con and the p-u mixing cases. So it may not be very surprising that the physics 
solution of any system corresponds to the one where the modula of the amplitudes take 
the minimal values among all the solutions. We call this selection rule "the minimal 
amplitude conjecture". 

Such a rule can be comprehended readily from a simple "economic" principle. One 
would expect that a solution in which the physics observable is produced via two or 
more very large amplitudes with strong destructive interference is less "economic" than 
the one via small amplitudes but with constructive interference. In the latter case, we 
always find maximum constructive interference among amplitudes. With this economic 
"minimal amplitude conjecture", we can determine the unique physics solutions of all the 
aforementioned examples. This will make the life much simpler. 

7 Conclusion 

In this Letter, we put forth two important problems: multiple solutions in fitting 
experimental data and the selection of the physics solution among many possible solutions. 
As a matter of fact, more and more cases with multiple solutions appear in recent data 
analyses and even more may have been overlooked previously. Here, we strongly suggest 
that all possible solutions be found out and reported in the future analyses, in order to 
avoid a babel of arguments and misleading theoretical deductions. 

When confronting with multiple solutions, one has to make a choice. We propose a 
"minimal amplitude conjecture": the solution with minimal modula of the amplitudes 
is the physics one. Such a rule ensures a unique solution is singled out from multiple 
solutions. This conjecture could be treated as a principle which should be tested experi- 
mentally by the future data analyses. 
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